SECOND ORDER STABILITY FOR THE MONGE- AMPERE EQUATION 
AND STRONG SOBOLEV CONVERGENCE OF OPTIMAL TRANSPORT MAPS 
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Abstract. The aim of this note is to show that Alexandrov solutions of the Monge- Ampere equation, 
with right hand side bounded away from zero and infinity, converge strongly in W^^^ if their right hand 
side converge strongly in I/ioc- As a corollary we deduce strong W^^^ stability of optimal transport 
maps. 



1. Introduction 

Let Q C be a bounded convex domain. Recently [3] the authors showed that convex Alexandrov 
solutions of 

( det D'^u = f inn 

(1.1) <^ 

I li = on dQ, 

with < A < / < A, are W^^^{fl). More precisely, they were able to prove uniform interior LlogL- 
estimates for D'^u. This result has also been improved in [H [8], where it is actually shown that 
u G Wf^;!in) for some 7 = 7(n, A, A) > 1: more precisely, for any $7' CC Q, 

(1.2) / \D^u\^ < C{n, A, A, n'). 

Jn' 

A question which naturally arises in view of the previous results is the following: choose a sequence 
of functions fk with A < < A which converge to / strongly in L^^^^Q), and denoted by Uk and u 
the solutions of (jl.lh corresponding to fk and / respectively. By the convexity of Uk and u, and the 
uniqueness of solutions to (jl.l|) . it is immediate to deduce that Uk u uniformly, and Vufc 
in Lj'j^^(r2) for any p < oo. However, what can be said about the strong convergence of D'^Uk"! Due 
to the highly nonlinear character of the Monge-Ampere equation, this question is nontrivial. (Note 

2 1 

that weak W^^^ convergence is immediate by compactness, even under the weaker assumption that fk 
converge to / weakly in L\^^{Q).) 

The aim of this short note is to prove that actually strong convergence holds. In fact our main 
result is the following: 

Theorem 1.1. Let Qk C M" be a family of convex domains, and : 11^ ^ R &e convex Alexandrov 
solutions of 



(1.3) 



det D'^Uk = fk in Qk 
Uk = on dQ,k 



2 
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with < A < /fc < A. Assume that converge to some convex domain Q in the Hausdorff distance, 
(ind fkXfik converge to f in Ll^^{Q). Then, if u denotes the unique Alexandrov solution of 

( det D'^u = f inn 
1 u = on dQ, 

for any Q' C ^ we have 

(1.4) — ^^||^y2,l(Q') — > as k ^ oo. 

(Obviously, since the functions Uk are uniformly hounded in W'^''^{Q.'), this gives strong convergence 
in W'^''^ (fi') for any 7' < 7. J 

As a consequence of the previous theorem we can prove the fohowing stabihty result for optimal 
transport maps: 

Theorem 1.2. Let Qi,Q2 C M" be two bounded domains with Q2 convex, and let fk,gk be a family of 
probability densities such that < A < fk,9k ^ ^ inside and il.2 respectively. Assume that fk^f 
in L^{Qi) and gk ^ g in L^{^l2), and let : f^i — )• Q2 (resp. T : f^i — )• Q2) be the (unique) optimal 
transport map for the quadratic cost sending fk onto gk (resp. f onto g). Then Tk ^ T in Wj^'J {^i) 
for some 7' > 1. 

We point out that, in order to prove (jl.4p and the local W^'^ stability of optimal transport maps, the 
interior L log L-estimates from [3] are sufficient. Indeed, the VF^''''-estimates are used just to improve 
the convergence from W^^^ to VFj^'J with 7' < 7. 

The paper is organized as follows: in the next section we collect some notation and preliminary 
results. Then in Section [3] we prove Theorem I l.H and in Section H] we prove Theorem 11.21 

Acknowledgments: We thank Luigi Ambrosio for stimulating our interest on this problem. We 
also acknowledge the hospitality at the University of Nice during the workshop "Geometry meets 
transport", where part of the present work has been done. AF has been partially supported by NSF 
Grant DMS-0969962. Both authors acknowledge the support of the ERG ADG Grant GeMeThNES. 

2. Notation and preliminaries 
Given a convex function li : — )> M, we define its Monge-Ampere measure as 

^iu{E) := \du{E)\ V^cJ^Borel 

(see [71 Theorem 1.1.13]), where 

du{E) := U du{x). 

Here du{x) is the subdifferential of u at x, and |F| denotes the Lebesgue measure of a set F. In case 
u G Cj^'c, by the Area Formula [5l Paragraph 3.3] the following representation holds: 

= det D^u dx. 

The main property of Monge-Ampere measure we are going to use is the following (see (71 Lemma 
1.2.2 and Lemma 1.2.3]): 
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Proposition 2.1. Let Uk : Q ^ M be a sequence of convex functions converging locally uniformly 
to u. Then the associated Monge-Ampere measures fiu^ converges to fiu in duality with the space of 
continuous functions compactly supported in Q. In particular 

k—^oo 

for any open set A C il. 

Given a Radon measure z/ on M" and a bounded convex domain Q C M*^, we say that a convex 
function u : — ?> M is an Alexandrov solution of the Monge-Ampere equation 

det D^u = V in 

if ^u{E) = i^{E) for every Borel set E C ^l. 

If v : i7 ^ M is a continuous function, we define its convex envelope inside as 

(2.1) r„(x) := sup{£{x) : e<vmn, i affine}. 
In case is a convex domain and v G C^(r2), it is easily seen that 

(2.2) D'^v{x) > for every x e {v = r^}nn 

in the sense of symmetric matrices. Moreover the following inequality between measures holds in Q: 

(2.3) fir, < det i:>^i;l{„=r„} dx, 

(here 1^; is the characteristic function of a set E). To see this, let us first recall that by [71 Lemma 
6.6.2], if xq G \ {Ft, = v} and a G dTi,{xQ) then the convex set 

{x G S7 : r^(x) = a ■ {x - xo) + r^,(xo)} 

is nonempty and contains more than one point. In particular 

dTy{n \ {r„ = v}) C {p G M" : there exist x,y e n, x ^ y and p e dTy{x) n dr^{y)}. 

This last set is contained in the set of nondifferentiability of the convex conjugate of u, so it has zero 
Lebesgue measure (see [3 Lemma 1.1.12]). Hence 

(2.4) |ar„(J7\{r„ = t>})| =0. 

Moreover, since u G C^{il), for any x G {Ty = v} (iQ it holds dT^(x) = {'Vv{x)}. Thus, using (12. 4p 
and (|2.2p . for any open set j4 CC we have 

//r„ {A) = I dTy {A n {r„ =v})\ = \Vv{An {F, = v}) I 

< / \detD^v\ = / detD'^v. 

JAn{ry=v} JAn{ry=v} 

(The inequality above follows from the Area Formula in [5} Paragraph 3.3.2] applied to the map 
S/v.) This proves (|2.3p . 

We recall that a continuous function v is said to be twice differentiable at x if there exists a (unique) 
vector Vu(x) and a (unique) symmetric matrix V^ti(x) such that 

v{y) = v{x) + Vu{x) ■ {y — x) + -V^n(x)[y — x,y — x] + o{\y — xp). 

In case v is twice differentiable at some point xq G = F„}, then it is immediate to check that 

(2.5) V\(xo) > 0. 
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By Alexandrov Theorem, any convex function is twice differentiable almost everywhere (see for in- 
stance [H Paragraph 6.4]). In particular (I2.5P holds almost everywhere on {v = T^}, whenever u is a 
finite linear combination of convex functions. 

2 1 

Finally we recall that, in case v € W^^^ , then the pointwise Hessian of v coincides almost everywhere 
with its distributional Hessian (5l Sections 6.3 and 6.4]. Since in the sequel we are going to deal with 

2 1 9 

VFjq'j, convex functions, we will use D u to denote both the pointwise and the distributional Hessian. 

3. Proof of Theorem 11.11 

We are going to use the following result: 

Lemma 3.1. Let $7 C M" be a bounded convex domain, and u,v : ^ M. be two strictly convex 
functions such that fi^ = f dx and fj,y = gdx, with f,g£ Lj^^{Q,). Then 

(3.1) /^r„_„ < - 5'/")"l{n-.=r„_4 dx. 

Proof. In case u,v are of class C^{Q,), by (12. 2p we have 

< D'^u{x) — D'^v{x) for every x ^ {u — v = T^-v}, 

so using the monotonicity and the concavity of the function det^^" on the cone of non- negative sym- 
metric matrices we get 

Q<dei{D^u- D^v) < ((detD^^)^/" - (detD^t,)^/")" on {u - v = Tu-v}, 



which combined with (|2.3p gives the desired result. 

Now, for the general case, we consider a sequence of smooth uniformly convex domains increasing 
to n, two sequences of smooth functions and gk converging respectively to / and g in Li^^{Q,), and 
we solve 

j det D'^Uk = fk in ftk \ det D^Vk = gk in ftk 

]^Uk = u*pk on d^k, \vk = v*pk on (90fc, 

where pk is a smooth sequence of convolution kernels. In this way both Uk and Vk are smooth on Qk 
[Bl Theorem 17.23], and they converge locally uniformly to u and v respectively. Hence, also Tui^-v^. 
converges locally uniformly to r„_^. Moreover, it follows easily from the definition of contact set that 

(3-2) limsupl|„^_^^=r„,_,j < Mu-v=r^.,}- 

fc— >oo 

We now observe that the previous step applied to Uk and Vk gives 



n 



m < \^{det D^Uk)^^"' - {detD^VkY^'' j l{u^-v^=r^^.^^} dx, 

Thus, letting k —?■ oo and taking in account Proposition 12.11 and (j3.2p . we obtain (13. ip . □ 

Proof of Theorem \l.l[ The L^^^ convergence of Uk (resp. Vu^) to u (resp. Vu) is easy and standard, 
so we focus on the convergence of the second derivatives. 

Without loss of generality we can assume that is convex, and that fi' CC (since 0^ ^ O in 
the Hausdorff distance, this is always true for k sufficiently large). Fix e € (0,1), let r„„(]^_j)„^ be 
the convex envelope of n — (1 — e)uk inside (see (12. 1|) ). and define 

Al := {xen' : u{x) - (1 - e)uk{x) = r„_(i_^)„^ (x)}. 



STABILITY FOR THE MONGE-AMPERE EQUATION AND STRONG CONVERGENCE OF OPTIMAL MAPS 5 



Since Uk ^ u locally uniformly, r„_(i_£)„^ converges uniformly to T^u = (as u is convex) inside fi'. 
Hence, by applying Proposition 12.11 and (|3.ip to u and (1 — e)uk inside Q', we get that 



Jn' 

< liminf/ir„_(i_,) (!^') 

< liminf / f/^/" - (1 -e)/; 



We now observe that, since fk converges to / in Ll^^{Q,), we have 



/ (/V- - (1 - e)fl 
Jn'nAi ^ 











' Jn'nAi 





as A; — )> cxD. Hence, combining the two estimates above, we immediately get 



/ < liminf 



f, 



or equivalently 



lim sup / / = 0. 

k^oo Jn'\Al 



Since / > A inside ft (as a consequence of the fact that /fc > A inside 17^), this gives 
(3.3) lim \n'\Al\ =0 Ve G (0,1). 

k—^oo 

We now recall that, thanks to [B O [H [8], both u and (1 — e)uk are strictly convex and belong to 
W'^'^{Q'). Hence we can apply (j2.5p to deduce that 

D^u-{l-e)D'^Uk>0 a.e. on^|. 

In particular, by (|3.3p . 

\{D^u<{l-e)D^Uk}nn'\^0 as/c^oo. 
By a similar argument (exchanging the roles of u and Uk) 

\{{l-e)D'^u>D'^Uk}nn'\->0 as/c^oo. 

Hence, if we cah := |x G f]' : (1 - e)D'^Uk < D^u < j^D'^Uk^ it holds 



lim \n'\Bl\ = Ve G (0,1) 

fe— >oo 
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Moveover, by (ll.2p applied to both Uk and u, we 

\D^u-D^Uk\= / \D'^u- D^Uk\+ / \D'^u-D'^Uk\ 



n' Jn'nBf, Jn'\Bf, 



[ \DW + \\D\-D\k\\LHn')\^'\Bl\'~'^^ 
- ^ Jn' 

<C(^ + If]'\i?ir-V7). 



Hence, letting first A; — )■ oo and then sending e — )• 0, we obtain the desired result. □ 

4. Proof of Theorem 11.21 

In order to prove Theorem 11.21 we will need the following lemma (note that for the next result we 
do not need to assume the convexity of the target domain) : 

Lemma 4.1. Let Qi,Q2 C M" be two bounded domains, and let fk,gk be a family of probability 
densities such that < A < /fe, (/^ < A inside Qi and Q2 respectively. Assume that fk^fin L^{Qi) 
and Qk ^ g in L^{yt,2), and let : — )• ^2 (resp. T : ^ Q2) be the (unique) optimal transport 
map for the quadratic cost sending fk onto gk (resp. f onto g). Then 

• ^ inL\ni). 



QkoTk goT 

Proof. By stability of optimal transport maps (see for instance 0| Corollary 5.23]) and the fact that 
/fc > A (and so / > A), we know that — ?• T in measure (with respect to Lebesgue) inside 

We claim that g oT^ — t- (7 o T in L^{Qi). Indeed this is obvious if g is uniformly continous (by the 
convergence in measure of Tfc to T) . In the general case we choose G C(r22) such that H^— 5'r;||Li(n2) — 
rj and we observe that (recall that fk,f>^,9ki9^-^y and that by definition of transport maps we 
have r#/fc = gk, T#f = g) 

[ \goTk-goT\<[ \gr,oTk-gnoT\+ [ \gr, oTk - g oTk\^ + [ \g^oT-goT\^ 

jQi JQi JQi ^ JUi 

= / \9v°Tk- gr,oT\+ \9v-9\^ + 

f , , ^ 

</ \9v°Tk- 9v°T\ + '2-i]. 



\9v-9\j 



Thus 



lim sup / \g o Tk — g o T\ < 2— rj, 
Jfii 



fc-s-oo -/o, A 
and the claim follows by the arbitrariness of rj 



If instead of p.2|) we only had uniform LlogL a-priori estimates, in place of Holder inequality we would apply the 
elementary inequality t < St log{2 + t) + e^^^ with t = \D'^u — D'^uu\ inside fi' \ and we would let first fc — > 00 and 
then send 5, £ ^ 0. 
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Since 



gk°Tk-goT\< I \gk oTk - g oTk\^ + [ \goTk-goT\ 

\gk-g\^+ \goTk-goT\ 
n2 -'Hi 



< \\\gk - g\\LHn2) + [ \g°Tk-goT\ 



from the claim above we immediately deduce that also gk o ^ g o T in L^{Q,i). 
Finally, since gk,g > A and / < A, 



fk 



f 



gk°Tk goT 



< 



fk - f 



gk ° Tk 



1 



1 



1 



5fc o Tfc goT 
gkoTk- go T\ 
gkoTkgoT 



<T\\fk-f\\mn,)+^ 

< -^Wfk - /llLi{f^i) + ^llS'fc °Tk - go ^llLi(ni)) 



from which the desired result follows. 



□ 



Proof of Theorem \1.2l Since are uniformly bounded in W^''^{^1'^) for any CC it suffices to 
prove that Tfc^T in VFi^^'^^(17i). 

Fix xq G Qi and r > such that Br^xo) C By compactness, it suffices to show that there is an 
open neighborhood Uxq of xq such that Uxq C Br{xo) and 



rl + iVTfc - VTI 0. 



It is well-known [2] that (resp. T) can be written as Vn^ (resp. Vit) for some strictly convex 
function Uk '■ Br{xQ) — >■ M (resp. u : Br{xQ) — >■ M). Moreveor, up to remove an additive constant 
(which will not change neither Tj. nor T), one may assume that Ufc(xo) = ^^(a^o)- 

Since Tk = Vuk are bounded (as they take values in the bounded set ^2)1 by classical stability of 
optimal maps (see for instance [U Corollary 5.23]) we get that Vuk — > in L\^^{Br{xQ)). (Actually, 
if one uses [2], Vuk are locally uniformly Holder maps, so they converge locally uniformly to Vn.) 
Hence, to conclude the proof we only need to prove the convergence of D'^Uk to D"^ u in a neighborhood 
of Xq. 

To this aim, we observe that, by strict convexity of n, we can find a linear function l{z) = a ■ z + b 
such that the open convex set Z := {z : u{z) < u{xo) + i{z)} is non-empty and compactly supported 
inside Bj./2{xo). Hence, by the uniform convergence of Uk to u (which follows from the Lj^^ convergence 
of the gradients, the convexity of Uk and u, and the fact that Uk{xo) = u{xo)), and the fact that 
is transversal to £ on dZ, we get that Zk := {z : Uk{z) < Uk{xo) + i{z)} are non-empty convex sets 
which converge in the Hausdorff distance to Z. 

Moreover, by [2] the maps Vk := Uk — i solve in the Alexandrov sense 



'det D'^Vk - 
Vk = 



in Zk 
on dZk 



8 
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(here we used that the Monge- Ampere measures associated to Vk and Uk are the same). Therefore, 
thanks to Lemma [4. II we can apply Theorem ll.il to deduce that D'^Uk — > D^u in any relatively compact 
subset of Z, which concludes the proof. □ 
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